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Abstract — This paper proposes novel joint channel estimation 
and heamforming approach for mnlticell wideband massive mnl- 
tiple input multiple output (MIMO) systems. Using our channel 
estimation and beamforming approach, we determine the number 
of cells Nc that can utilize the same time and frequency resource 
while mitigating the effect of pilot contamination. The proposed 
approach exploits the multipath characteristics of wideband 
channels. Specifically, when the channel has L multipath taps, 
it is shown that Nc < L cells can reliably estimate the channels 
of their user equipments (UEs) and perform beamforming while 
mitigating the effect of pilot contamination. For example, in a 
long term evolution (LTE) channel environment having delay 
spread Td = 4.69/i second and channel bandwidth B = 2.5MHz, 
we have found that L = 18 cells can use this band. In practice, Td 
is constant for a particular environment and carrier frequency, 
and hence L increases as the bandwidth Increases. The proposed 
channel estimation and beamforming design is linear, simple to 
implement and significantly outperforms the existing designs, and 
is validated by extensive simulations. 

Index Terms — Beamforming, Channel estimation. Massive 
MIMO, Pilot contamination, Rayleigh Quotient 

I. Introduction 

Massive multiple input and multiple output (MIMO) is one 
of the promising techniques to satisfy the gigabit connectivity 
requirements of the future 5G network [1]. In a multiuser and 
multicell setup, massive MIMO can be deployed both at the 
base stations (BSs) and user equipments (UEs). However, since 
the energy consumption and cost of transceiver devices in¬ 
crease as the number of antennas increases, deploying massive 
antennas at the UEs is usually infeasible from practical aspect 
particularly at microwave frequency bands. This motivates the 
current paper to consider a massive MIMO system where each 
BS (UE) has massive antennas (single antenna) as in [2], [3]. 

The potentials of a massive MIMO system can be exploited 
by performing beamforming which depends on the availability 
of accurate channel state information (CSI). However, in prac¬ 
tice, the channel between each BS and UE is estimated from 
orthogonal pilot sequences which are limited by the coherence 
time of the channel [3]-[5]. Eor a massive MIMO system, 
one can possibly estimate and exploit the maximum number 
of channel coefficients using time division duplex (TDD) 
approach by sending pilots from UEs [5], [6]. By doing so, 
each BS will estimate, and utilize the channel coefficients both 
for the downlink and uplink channel data transmissions. In a 
multicell setup where there is no coordination between cells. 


each BS determines its number of UEs from the coherence 
time of the channel only and the UEs in each cell utilizes the 
same pilot sequences. The reuse of pilot sequences of several 
co-channel cells will create a so called pilot contamination 
which is shown to be a fundamental performance bottleneck 
for multicell massive MIMO systems [l]-[3], [5]-[9]. 

A number of approaches are proposed to address pilot con¬ 
tamination. Eigenvalue decomposition (EVD) based channel 
estimation approach is proposed in [7] where it is demon¬ 
strated that this method achieves better estimation accuracy 
than that of the conventional least square (LS) and minimum 
mean square error (MMSE) estimators. A coordinated channel 
estimation approach is proposed in [2] where it is shown that 
pilot contamination can be vanished when UEs have special 
channel covariance matrices. In [8], pilot contamination elim¬ 
ination approach is proposed by allowing pilot transmissions 
both in the uplink and downlink channels. Successive pilot 
transmission approach has been proposed in [9] to eliminate 
pilot contamination. This paper utilizes consecutive pilot trans¬ 
mission phases in which each BS stays idle at one phase and 
repeatedly transmits pilot sequences in other phases. 

In [4], [10], time shifted channel estimation approach for 
orthogonal frequency division multiplexing (OEDM) based 
frequency selective channels is proposed. In such approach, 
first the cells are grouped in some predefined way. Then, only 
one group of cells send pilot signal sequentially while all 
the other groups of cells transmit their downlink data signal. 
Einally, all cells transmit their uplink data on the remaining 
time slot. In this approach, however, since the number of cells 
in each group could be more than one, large coherence time 
is required to cancel pilot contamination in each group of 
cells when every cell serves the same number of UEs like 
in the conventional LS and MMSE approaches (for example 
like in [3], [4]). Erom these discussions, one can notice that 
the problem of pilot contamination has been considered by 
several papers. However, some of the aforementioned papers 
address pilot contamination by exploiting particular channel 
covariance information which may not be hold in all scenarios, 
and the rest of them utilize extra pilots compared to [3], [4]. 

The current paper examines a pilot contamination problem 
for multicell massive MIMO systems with frequency selective 
channels. Eor a given pilot duration, we assume that the 
number of UEs served by each BS is the same as that of 
[3], [4]. Eurthermore, the channel covariance matrix of each 


UE is assumed to have ’’arbitrary structure”. Under these 
settings and assumptions, we propose a novel joint channel 
estimation and beamforming approach. The proposed approach 
exploits the inherent multipath behavior of frequency selective 
channels. Specifically, when the channel has L multipath 
taps, it is shown that Nc < L cells can reliably estimate 
the channels of their UEs and perform beamforming while 
efficiently mitigating the effect of pilot contamination. 

In other words, when < L, it is shown that the signal to 
interference plus raise ratio (SINR) of each UE’s sub-carrier 
increases as the number of antennas at each BS N increases 
(i.e., the SINR grows indefinitely as N ^ oo). In an exemplify 
long term evolution (LTE) channel environment having delay 
spread = 4.69p seconds and channel bandwidth B = 
2.5MHz, L = 18 cells can use this frequency band while 
efficiently mitigating the effect of pilot contamination. In 
practice, is constant for a particular environment and carrier 
frequency, and hence L increases as the bandwidth increases 
[II]. The proposed design is linear, simple to implement and 
significantly outperforms the existing designs, and is validated 
by extensive simulations. 

This paper is organized as follows. Section II discusses the 
system and channel models. In Section III, a summary of pilot 
contamination and the objective of the paper is explained. In 
Sections IV - VI, the proposed algorithm and its performance 
analysis is presented briefly. In Section VII, simulation results 
are presented. Einally, Section VIII draws conclusions. 

Notations: In this paper, upper/lower-case boldface letters 
denote matrices/column vectors of appropriate size. Identity 
matrix, singular value decomposition, expectation and absolute 
value are denoted by I, SVD(.), E(.) and |.|, respectively. 
Greater than or approximately equal to is denoted by 


II. System and Channel Model 

We consider a multiuser and multicell system with cells 
(i.e., BSs) where each UE and BS are equipped with 1 and N 
antennas, respectively. We consider transmission scheme with 
symbol period Tg and channel environment with maximum 
delay spread T^'. The total transmission bandwidth becomes 
B = i?o(l+a), where Bq = ^ is the Nyquist bandwidth and 
a is the excess bandwidth [13], [14]^. Eor these settings, the 
number of multipath channel taps L is approximated as L = 
^ [11] (see page 143) and the channel coefficients between 
the fcth UE in jth cell to the ith BS nth antenna can be denoted 
as [12], [14] 


^kjin — [hkjinl^ ‘ ‘ ‘ 


( 1 ) 


where — y/9^jj^kjini with , Vfc, j, z, n, ( (^(]kji) 

fast fading (distance dependent slow fading) coefficients. 

In most systems having multipath channels, OEDM based 
transmission is adopted. Eor such a transmission, the channel 
coefficient of each sub-carrier has practical importance. The 


* The transmission environment has a maximum delay spread of Td which 
has typical values = {4.7, 5.2}/rs for urban cells and = 16.7/rs for 
very lai'ge multi-cell MIMO systems in the existing LTE network [3], [12]. 

-The exact value of a varies from one application to another (for example, 
0.2 < o < 0.35 is used in [13]). 


channel coefficient between the A:th UE in )th cell to the zth 
BS nth antenna in sub-carrier s can be computed from (1) as 
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kjins 
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kjin 


( 2 ) 

where = [1, 

with M as the fast Eourier transform (EET) size of the OEDM. 


III. Pilot Contamination and Objective 


The number of terminals served by a BS is proportional to 
the pilot duration Tp. Eor the given Tp, the channel coefficients 
of each UE can be estimated with or without OEDM pilot 
transmission. By denoting each OEDM duration as To and 
useful symbol duration as Tu, each BS can serve the following 
maximum number of UEs (see (4) of [3] and Section II of [4]): 


TpTu = 

Td To ~ Td L 


(3) 


where Np = 7^4 is the number of samples acquired in Tp. If 
each UE applies non OEDM pilot transmission, it is possible 
to serve K UEs [3]. 

If a BS serves K UEs (i.e., use Tp to estimate the channel 
coefficients of these UEs), we can understand from the work 
of [3] that only one cell can serve its UEs without experiencing 
any pilot contamination. In particular, when hkjini,yk, j, i, n, I 
are modeled as independent and identically distributed (i.i.d) 
zero mean circularly symmetric complex Gaussian (ZMCSCG) 
random variables each with unit variance (i.e., hkjini ~ 
CJ\f(0,1),yk, j,i,n,l). And if there are multiple cells where 
the BS of each cell has sufficiently large number of antennas 
(i.e., N oo), the effect of noise vanishes and the SINR of 
the /cth UE in cell i and sub-carrier s is given as 


SINRk^s = 
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(4) 


As we can see from this equation, SINRkis is bounded 
even if N ^ oo. However, one can understand from MIMO 
communication that increasing N should help to improve 
SINRkis if pilot contamination is mitigated or cancelled. This 
motivates us to consider the following objective: 

Eor the given B, Tp, To, Tg, Td and r„, the channel vector 
between any two UEs are uncorrelated and each BS serves K 
UEs (i.e., the same setting as in [3], [4]), how many cells can 
utilize the same time-frequency resource while ensuring that 
SINRkis increases as N increases (i.e., achieving unbounded 
SINRkis as N ^ oo which is the same as mitigating or 
possibly canceling the effect of pilot contamination)? 


IV. Proposed Channel Estimation and 
Beamforming 

This section discusses the proposed channel estimation and 
beamforming approach. In this regard, we assume that the 
pilots are transmitted without OEDM whereas, the data sym¬ 
bols are transmitted using OEDM scheme. The current paper 
considers that both pilot and data transmissions take place in 
the uplink channel. We also assume that TV), (which will be 
determined in the sequel) cells utilize the same time-frequency 
resource while mitigating the effect of pilot contamination. 
Here, we consider that the zth cell is the target cell. 




A. Channel Estimation 

This subsection provides the proposed channel estimation. 
For the considered settings, the received signal at the nth 
antenna of the ith BS can be expressed as 
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Xi = [Xii,X 2 i,- • • ,x^-] and Wi„ G is the additive 

noise vector during pilot transmission phase where its entries 
are assumed to be i.i.d ZMCSCG random variables each with 
variance and Xkj € qN^xl ^ Toeplitz matrix, i.e.. 
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with Xfej = [xkj^,Xkj 2 - ■ ■ the pilot symbols 

transmitted from the fcth UE in cell j which will be designed 
in Section V. As can be seen from (5), the ith BS experiences 
pilot contamination due to the term Xihun [3]. For an 
arbitrary pilot, one can also notice from this equation that 
Xihiin = 0 is ensured when Nc = 1 (i.e., like in [3], [4]). For 
such setting, one can apply MMSE or FS method to estimate 
hji„ from Tin- 

In the following, we provide the proposed channel estima¬ 
tion approach. In this regard, we introduce a vector Vkis G 
C^pXi J.Q estimate the sth sub-carrier channel of the fcth UE 
in cell i which will be designed in Section V. Using this linear 
combination vector, we express the estimate of hkuns as 

tv, _ 

hkiins — ^in^kis — (X^h^iyi -f ^ ^ X^jYlj^n “t” Stkis- 

(7) 

As will be clear in the next section, this estimated channel 
helps us to increase the number of cells more than one. 


B. Beamforming 

As discussed above, we consider the uplink channel for data 
transmission using OFDM approach. Since data transmission 
takes place using OFDM technique, the received signal of 
each sub-carrier can be obtained independently at each antenna 
of all BSs. To this end, the ith BS nth antenna receives the 
following uplink signal in sub-carrier s (Uins) 


K Af, 


Uins — ^kjinstlkjs “t" ‘^ins 

fc=li=l 


where Wms is the noise sample at the ith BS nth antenna 
and sub-carrier s during data transmission and dkjs is the data 
symbol transmitted in sub-carrier s of the fcth UE in cell j. The 
overall received signal at the ith BS sth sub-carrier becomes 

k Af, 

Yis ^ ^ ^ ^ ^kjist^kjs ^is 

k=lj=l 

k jv, 
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where yi^ = [yiis,yt2s,-■ ■ ,yiNsV' = 

[wiis,Wi 2 s, -'' It is assumed that each entry 

of Wis is assumed to be i.i.d ZMCSCG random variable with 
variance Now, let us assume that we are interested to 
estimate dkis using a beamforming vector akus G as 

dkis = ak^isy^s ( 8 ) 

k Af, 

— ^kiisO^kiistikis T ^ ^ ^ ^ ^mjisdmjs ^is)- 

m=l j=l,(m,j)^(k,i) 


Using this beamformer, dkis will have the following SINK 
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As can be seen from this equation, the achievable fkis depends 
on akiis- Furthermore, one can utilize different beamforming 
methods to design akus- In a massive MIMO system, simple 
beamforming methods such as maximum ratio combining 
(MRC) and zero forcing (ZF) are close to optimal which 
motivates us to choose akus to be the MRC beamformer [3], 
[5], [ 6 R That is akus = hkiis [^kiilst ^kii2s ‘ ‘ ‘ i ^kiiNs] 5 
where hkuns is the estimated channel defined in (7). With this 
beamforming vector, we will have 
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with = E{H;;,^-iHj„f,*fjH;„H^^-J and = 

E{H^^jf*fjH^jj} (see Appendix A of [15] for the details). 

The beamforming vector akus exploits the estimated chan¬ 
nel (7) which is a function of x^j, Vm, j and TV),. Hence, 
the achievable fkis depends on these variables that need to be 
optimized which is the focus of the next section. 


V. Optimization of 'tCmj and ^tkis 

This subsection discusses the optimization of x^j and 
Nc- In this regard, we examine two problems: In the first 
problem, we determine the number of cells Nc, and feasible 
Vkis and 'x.mj such that the multi-cell system is able to mitigate 
pilot contamination. In the second problem, we re-optimize 
Vfeig to further maximize fkis for fixed Nc and Xmj,Vm,j. 










A. Determination of Nc, and Feasible 'Vkis tind 

This subsection determines Nc, and feasible Vkis and 
such that the multicell system is able to mitigate the effect of 
pilot contamination. Specifically, we determine these variables 
while ensuring that fkis increases as N increases (i.e., fkis 
grows to infinity as N —t oo) in the following Theorem. 

Theorem 1: For arbitrary channel covariance matrices given 
in (10), fkis of (9) increases as N increases (i.e., the effect of 
pilot contamination is mitigated) when Nc < L, Xmj,Vm, j 
are selected such that s are full row rank matrices and 

Vfcis is designed from the solution of 

Ukiis (11) 

where = [quisj ft 21 is) ' ' ' ^KUs^ C[l 2 isi 9t22is) ' ' ' , Q.K 2 is’ 
iQ/fwis]’ ^mjis ~ ^raji^rny 

= SVD(q£), U,, e and 

Uis G (^NpxNj, unitary matrices, is the row 

vector of corresponding to q^jjg and G is 

a diagonal matrix containing the singular values of Q^. 

Proof: See Appendix B of [15]. ■ 

For arbitrary covariance matrices C() of (10), one approach 
of satisfying full rank is by selecting Xmj,Vm,j so that 
they will be uncorrelated to each other. Furthermore, in order 
to maintain full rank for all sub-carriers, x^^ , Vm, j may 
need to have a ’’white noise” like sequence. This motivates 
us to select Xn from noise like deterministic (or random) 
sequence, and then constmct Xmj,V(TO, j) f (1,1) by shifting 
the elements of Xn with N^ sized FFT matrix. With such 
selection, however, we have noticed that the average rate 
achieved by one sub-carrier differs from the other sub-carrier 
and this rate imbalance may not be desirable in practice. For 
this reason, we choose Xn randomly from a given set of 
sequences where these sequences are known a priori to all 
BSs and UEs. Upon doing so, we have observed that each 
sub-carrier of a UE achieves the same average rate. 

One can understand from the result of Theorem 1 that if 
Nc > L, the proposed channel estimation and beamforming 
approach will achieve ’’bounded fkis" even though N ^ oo 
(i.e., if Nc > L, we will experience pilot contamination as 
will be demonstrated in the simulation section). 

B. Re-optimization of stkis 

Erom the above subsection, we are able to determine Nc, 
and provide feasible Vkis and x^j while ensuring that the 
effect of pilot contamination is mitigated. Eor the given Nc, it 
is also possible to optimize Vkis and to further increase 
fkis- However, jointly optimizing Vkis and x^j to maximize 
fkis is still complicated. Eor this reason, we re-optimize Vkis 
only to maximize fkis for fixed Nc and x^j as follows. 

max fkis- (12) 

Vkis 

This problem is a Rayleigh quotient optimization problem 
where generalized eigenvalue solution approach can be applied 
to get the optimal Vkis [13], [16], [17]. 

The next issue is what is the relation between the solutions 
obtained from (11) and (12)? In the following, we address 
these issues; Denote the solution obtained by solving (12) as 


and its corresponding SINR as As this problem is 
a Rayleigh quotient, is a global optimal solution. Hence, 
the SINR obtained by any 'Vkis ^ can not be higher than 
that of which leads to (after some mathematical steps) 

iLs > ikrs (13) 


where 


N- 


'ykis — 


wherea = vL(Ef=iE:= 

La'^Ckiis)'Vkis, b = 

E 

(T^I)vfeis- As can be seen from this expression, a and b are 
constant terms which are not dependent on N. Hence, for 
fixed Nc, L and K, one can expect that the SINR obtained 
by utilizing Vkis of (11) and (12) will be almost the same for 
large N. However, the solutions of these two problems may 
have different SINRs when N is not sufficiently large. This 
fact has been verified in the simulation section. 

We would like to stress here that the current paper utilizes 
pilot transmission without OEDM as shown in (5). However, as 
clearly seen from this equation, performing pilot transmission 
without OEDM ’’alone” does not bring any new result to 
increase Nc- Nevertheless, by performing pilot transmission 
without OEDM, data transmission with OEDM, carefully 
selecting x^j, and introducing and optimizing the variable 
'Vkis, we are able to increase the number of cells to Nc = L. 
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VI. Performance Analysis 
In a practical wireless system, the rate achieved by dkis (i-e., 
Rkis) which is related to fkis has an interest. Eurthermore, 
Rkis will be maximum when each of the BSs has perfect CSI 
of its own UEs. This motivates us to examine Rkis of the 
proposed approach and that of perfect CSI scenario. In this 
regard, we consider the following theorem. 

Theorem 2 When hkji are modeled as in (4), tL » 
7 fels >> 1 Qkii ^ Praji (the case in practice), we have 


R, = RZs - RZs 


< l0g2(f?) = Co,= 


Rll 


RZs 


(15) 


N—¥oo 


where R^ and RH are the rate achieved by perfect CSI 
scenario and proposed approach, respectively and the approx¬ 
imation is because rj is independent of N and is given as 


Lb 


V = 
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Vkis— soln. of (11) 


with b (^^(m,j)^(fc, 2 ) [^^(ii,'u)^(m,_ 7 ) ^9mji9uvi^rnj-^rnj 

+ + 0-^(E(„j) + (r^l))vkis- 

Proof: See Appendix C of [15]. ■ 

We would like to mention here that the results of the current 
paper are provided for the MRC beamforming approach in 
the uplink channel only. However, the approach of this paper 
can also be extended for other beamforming techniques and 
downlink channel. Eor such cases, Nc and x,„j are the same 
as that of the current paper. However, the optimal solution of 







Vfcis may not be necessarily the same for all beamforming 
types and channels. The detailed extension is omitted due to 
space constraint but it can be found in [15]. 

Example; For better insight about the usefulness of the pro¬ 
posed method, we provide numerical examples by employing 
the parameters used in the LTE network. For this network, 
different possible bandwidths are utilized (see [12] for more 
details). Here we consider the case where B = 2.5MHz with 
Td = 4.69/rs, coherence time = 0.5ms, Tp = ^ and 
Tg = = 0.2605/j,s [3], [12]. With these settings, we get 

iT = ^ « 53, TV, = L = ^ « 18. (16) 

One can notice from (16) that our approach can serve LK « 
954 UEs while mitigating pilot contamination. This is similar 
to scaling the number of UEs using the same resource by 
L. Note that if some resources are used for computational 
purpose, K will be reduced slightly [3]. 


VII. Simulation Results 


This section presents simulation results by considering a 
multicell system with Nc cells where each cell serves K = A 
UEs. In this regard, we set L = 4, gku = 1, gk2i = 0.9, gksi = 
0.6, gkii = 0.7, Vfc, Np = 16, M = 64, N = 2^» and k = 
1 = 1 (i.e., the desired user and cell) with Nq as an integer. 
We examine the average achievable rate of the UE’s sub-carrier 
(i.e., Rkis)- The SNR is defined as 7 = ,Vfc, i, 

All results are plotted by averaging 10000 channel realizations. 


A. Comparison of the Proposed and Existing Approaches 

This subsection compares the performance of the proposed 
approach with those of the existing ones. In the proposed 
approach, first we choose xn from a randomly selected unit 
energy quadrature phase shift keying (QPSK) signal. Then, 
we compute the remaining pilots just by shifting xn appro¬ 
priately with equally spaced conventional (fractional) discrete 
Eourier transform (DET) vectors. Eor the existing approaches, 
we design Xfci,Vi, fc in two methods; the first method uses 
orthogonal pilots like in the proposed approach and the second 
method applies full pilot reuse (i.e., x^^ = xj.j,Vi, j. A:). 

Eor fair comparison, we have utilized ^ = 1 (to resolve the 
multiplicative factor ambiguity) and binary phase shift keying 
transmitted signal of size 50 for the EDV approach of [7], 
and a full pilot reuse for [9]. Eig. 1 shows comparison of the 
proposed and existing approaches when 7 = OdB, s = 20 
and N(, = L. Erom this figure, we can observe that utilizing 
orthogonal pilots slightly improve the performances of the 
MMSE and LS methods. Eurthermore, the approaches of [7] 
and [9] also behave similar to the MMSE and LS approaches. 
However, the rate achieved by the proposed approach in¬ 
creases progressively with N. And, our approach achieves 
significantly higher Rkis compared to those of the existing 
approaches especially in a massive MIMO regime which is 
desirable. 


B. Validation of Theoretical Results 

This subsection validates the theoretical results provided in 
Sections V and VI. According to Theorem 1, if W, > 4, 
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Fig. i. Comparison of the rates achieved by the proposed and existing 
approaches. 


- A - Random (s=15, with (11)) 

0 Random (s=15, with (12)) 

A 

- - Random (s=27, with (11)) 

_B_ Random (s=27, with (12)) 


- A - Fixed (s=15, v . with (11)) 

^ * /f ** - 

' 'jp'/ 

- 4 . Fixed (s=15, with (12)) 

■ 

- Fixed (s=27, with (11)) 


- B - Fixed (s=27, with (12)) 



'A 


✓ 


. - 




2 4 6 8 10 12 14 16 

Normalized number of antennas (N^) 


Fig. 2. Comparison of the rates achieved by fixed and random pilots. 


increasing N indefinitely does not help to improve fkis- Also, 
the rate achieved by utilizing the solutions of (11) and (12) 
will be close to each other when N is very large and cq of 
(15) is derived for N ^ oo. These motivate us to perform 
simulation for this subsection by considering very large N. 

1) Effect of Pilots: This simulation examines the effect of 
pilots on the performance of the proposed algorithm. The pilot 
samples are designed by two approaches: The first approach 
utilizes Xii = + sZ—lh-f) where [zi Z 2 ] are the 

last 2Np samples generated from the ’’Maximum Length 
Sequence” ±1 bits of size 64 (i.e., fixed pilot assignment ap¬ 
proach). The second approach designs the pilots like in Section 
VII-A (i.e., randomly selected pilot assignment approach). We 
examine the rate Rkis (i-e., rate per user per sub-carrier) when 
7 = OdB and s = {15, 27} as shown in Eig. 2. As can be seen 
from this figure, in both sub-carriers, Rkis increases with N 
by utilizing either fixed or randomly generated pilots. 

On the other hand, the average Rkis may not be necessarily 
the same for all sub-carriers when fixed pilot assignment 
approach is utilized. However, the rate of these sub-carriers are 
the same when pilots are selected randomly which is expected 
(i.e., for given design approach of Vkis)- This simulation also 
exploits the fact that the rates achieved by employing the 
solutions of (11) and (12) are closer to each other when N 
is very large. However, for small to medium N, higher Rkis 

















VIII. Conclusions 



Fig. 3. The rates achieved by the proposed and perfect CSI scenarios. 



Fig. 4. Effect of Nc for the proposed approach. 


is obtained when Vkis is optimized using ( 12 ) which fits with 
the theory. In the following simulation, we utilize the v^is 
designed by ( 12 ) with random pilot assignment approach. 

2) Verification of cq: This simulation verifies the bound 
derived in (15) while utilizing 7 = OdB and s = 20 as 
shown in Fig. 3. From this figure, one can observe that the gap 
between the rate achieved by the proposed and that of perfect 
CSI scenario Rg becomes constant {Rg = cq « 6 b/s/hz) 
after N « 2^® which fits with the theory. On the other hand, 
Rg decreases as N decreases which is desirable in practice. 
For instance, if we deploy TV = 16 antennas, the rate loss 
is Rg < Ib/s/hz. We would like to mention here that, in a 
massive MIMO regime, we have observed almost the same 
behavior as in this figure for 7 = 20dB. This is simply because 
the effect of noise vanishes for sufficiently large N [3], [9]. 

3) Effect of Nc-' As discussed in Section V, the proposed 
design experiences pilot contamination if Nf. > L = A. This 
simulation validates this claim. To this end, we set s = 32, 
Nc = {5, 6 }, the channel gains of the UEs corresponding to 
the first four cells are the same as the first paragraph of this 
section, and the UEs of the fifth and sixth cells have channel 
gains gk 5 i = 0.6,gk6i = 0.75, Vfc. Eig. 4 shows the Rkis 
achieved for these settings. As can be seen from this figure 
and Eig. 3, Rkis drops quickly as we increase Nc from 4 to 
{5,6}. Eurthermore, in both Nc = {5,6} cases, Rkis is not 
increasing further after approximately N = 2^® (i.e., bounded 
rate) which is inline with the theory provided in Theorem 1. 


This paper proposes novel joint channel estimation and 
beamforming approach for multicell wideband massive MIMO 
systems. Using the proposed approach, we determine Nc uti¬ 
lizing the same time-frequency resource while mitigating the 
effect of pilot contamination for arbitrary channel covariance 
matrix of each UE. Specifically, when the channel has L mul¬ 
tipath taps, it is shown that Nc < L cells can reliably estimate 
the channels of their UEs and perform beamforming while 
mitigating the effect of pilot contamination. Eurthermore, 
whenever A} > L, the rate of each UE will be bounded even if 
N ^ 00 . These facts have been demonstrated using extensive 
computer simulations. The proposed joint channel estimation 
and beamforming approach is linear, simple to implement and 
significantly outperforms the existing approaches. 
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